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SUMMARY 


Stresses were calculated for finite-size orthotropic laminates loaded by 
a frictionless steel pin in a circular hole of the same diameter. These cal- 
culations were based on finite-element analyses for six laminates: quasi- 
isotropic [ 0°/±45O/90O] g, 0°' 90°, [0O/90°]g, [±45°]s, and [0°/±45°]g. Stress 
concentration factors, based on nominal bearing stress, were determined for wide 
ranges of the ratios of width to diameter w/d and edge distance to diameter 
e/d. For reference purposes, the infinite-laminate case was also analyzed for 
each laminate. 

Orthotropy had a significant influence on both the magnitude and location 
of the maximum tensile stress concentration on the boundary of the hole. For 
example, the stress-concentration factors for the infinite-laminate cases ranged 
from 0.82 to 2.16, compared with 0.87 for the quasi- isotropic laminate. The 
laminates with 0^ plies developed the peak tensile stress near the ends of the 
pin/hole contact arc. But the 90^ and ±45^ laminates had peaks where ply 
fibers were tangent to the hole. 

The finite widths and edge distances strongly influenced the tensile 
stress concentration. For the practical range w/d ^ 3, the peak tensile 
stresses were as much as 80 percent larger than the inf inite- laminate reference 
value. Even for rather wide laminates with w/d = 10, some of the peak tensile 
stresses were amplified by nearly 20 percent. For e/d ^ 3, these stresses 
were amplified by as much as 50 percent. In contrast, the finite width and 
edge distance had little effect on shear-out and bearing stress concentrations. 
For w/d Z 3 and e/d ^ 3, the peak shear-out stress was amplified by less 
than 25 percent and peak bearing stress by only 10 percent. 


INTRODUCTION 

The stress concentrations at fastener holes are widely recognized as having 
a degrading effect on structural fatigue life. In conventional metallic struc- 
tures, fatigue cracks commonly initiate at these stress-concentration sites. 

In recognition of this, the effects of stress concentrations on crack initiation 
and growth have been incorporated into current methods for predicting fatigue 
life. However, in composite structures, the stress concentration effects are 
much more complicated and are not well understood. For example, holes can pro- 
duce either composite delamination, matrix crazing, fiber failures, or various 
combinations of these. Furthermore, for a pin- loaded hole, these different 
types of damage can produce three different modes of failure: net-tension, 

shear-out, or bearing. Before analytical models can be developed for these 
various failures, the stress concentrations at loaded holes must be better 
understood. To date, the few stress analyses that have dealt with this problem 
have been limited to rather specific cases. The purpose of the present paper 
is to present some of the needed background on stresses at loaded holes in 
compos i t e laminates. 



This study was based on a two-dimensional finite-element analysis of 
stresses in an orthotropic laminate that is loaded by a frictionless steel pin 
in a circular hole of the same diameter. Six graphite/epoxy laminates were 
analyzed; quasi- isotropic [0°/±45°/90°]g, 0®, 90®, [0®/90®]g, [±45®]g, and 
[0®/±45®]s« Finite-element models were developed for wide ranges of width-to- 
diameter ratio w/d and edge-distance-to-diameter ratio e/d. First, each 
laminate was analyzed for different values of w/d using one large value of 
e/d. Next, a range of e/d was used with one large w/d value. 

The results ate presented in two ways: first as stress distributions near 
the hole and then as stress-concentration factors for ranges of w/d and e/d. 
Results for the orthotropic laminates were compared with the quasi-isotropic 
laminate which served as a reference case. Similarly, the w/d and e/d 
effects were evaluated relative to an infinite-laminate reference case. 
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SYMBOLS 

diameter of hole/ m 

Young's modulus in x~ and y-directions, respectively/ MPa 
edge distance/ m 
shear modulus/ MPa 

shear stress-concentration factor/ (^xy^max/^b 
Kgi^ for an infinite laminate 

tensile stress-concentration factor/ (°'09)max/Sb 
for an infinite laminate 

gross-section stress-concentration factor / (^60)max/Sg 

pin load/ N 

radius of hole/ m 

polar coordinates/ m and deg 

nominal bearing stress/ P/td/ MPa 

nominal gross-section stress/ MPa 

laminate thickness/ m 

width of laminate/ m 

Cartesian coordinates/ m 
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(^00) max 
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xy 


(^xy) max 


one-half the angle of contact between pin and hole, deg 
angular location for (c^00)niax hole boundary, deg 
Poisson's ratio 
stress, MPa 

polar stress components, MPa 

maximum tangential stress on hole boundary, MPa 
Cartesian shear stress, MPa 

maximum T^y on shear-out plane (y = R) , MPa 


FINITE- ELEMENT ANALYSIS 

Figure 1 shows the general problem analyzed in this study. The laminate 
had a frictionless hole which was loaded by a steel pin of the same diameter. 

The laminate reacted to the pin load by a uniform stress on its left end. The 
hole diameter was unity and the various combinations of w/d and e/d were 
generated by varying w and e in the finite-element models. The w/d 
effects were analyzed using a single, large value of e/d (e/d = 10). Simi- 
larly, the e/d effects were studied using w/d = 20. For all cases, the hole 
was at a distance of lOd from the left (loaded) end of the laminate. 

Table I presents the elastic constants for the six laminates analyzed. 

These laminates were selected to cover a broad range of composite properties. 

The 0° laminate properties (ref. 1) were used with laminate theory to calculate 
the elastic constants for the other laminates. Although the [ 0^/145^/90^] s 
laminate was called a quasi- isotropic laminate, its calculated elastic constants 
satisfy all the criteria for isotropy. Consequently, the results for the quasi- 
isotropic case were compared with isotropic solutions from the literature as 
a check on the adequacy of the finite-element analysis. 

The finite-element model for one-half of a 20d by 20d square laminate is 
shown in figure 2. To illustrate the model details, it is shown in four parts. 
This was the largest model analyzed (w/d = 20 and e/d = 10). Smaller values 
of w/d and e/d were obtained by removing elements from the width and 
unloaded end of the model, respectively. This model employed two types of 
finite elements. Linear-strain (6-node) elements comprised the 2.5d by 5d 
rectangular region around the hole. Uniform-strain (3-node) elements comprised 
the part of the model outside this region. The model was analyzed using the 
BEND finite-element computer program described in reference 2. 

Note that the pin is also modeled in figure 2. The pin was loaded at 
its center and was connected to the laminate by short, stiff spring elements 
(0.002d long) between adjacent nodes across the pin/laminate interface. Because 
these springs had no transverse stiffness, they transferred only radial loads 
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and thereby produced the desired frictionless interface. Also, these springs 
were used to allow the pin and laminate to separate along part of their inter- 
face. This was accomplished by assigning a zero stiffness to springs along the 
separation arc. References 3 and 4 show that this separation was independent 
of loading level for the special case of an initially snug fit (zero clearance) 
between the pin and laminate. For this reason, a snug fit was used in the 
present study. For each combination of material and configuration, the separa- 
tion angle was estimated and a computer run was made. When a spring was com- 
puted to have a tensile force, its stiffness was set equal to zero and that 
case was rerun. Such reruns were needed for only a few cases. 


RESULTS AND DISCUSSION 
Evaluation of Finite-Element Procedure 

The finite-element procedure was evaluated by comparing results for the 
quasi- isotropic laminate with results in the literature for an isotropic, 
infinite sheet containing a rigid pin with a snug fit. The appendix describes 
the present finite-element analysis used to approximate this infinite-sheet case. 
The results for this special case are presented in figure 3 as stress distri- 
butions around the hole. The stresses are normalized by the nominal bearing 
stress Sj 3 , which is defined as P/td, where P is the pin load, t is the 
laminate thickness, and d is the hole diameter. 

Figure 3 also includes a curve for Orr = -4 /tt cos 0 because several 
authors have assumed this cosine distribution for bearing stress, either 
directly or indirectly, by assuming a rigid-body displacement for the pin. The 
present results (solid curve) and the results of reference 5 (dashed curve) are 
in general agreement with the cosine distribution, except for the angle of 
contact between the pin and the laminate. The 90^ contact angle 9^ for the 
cosine distribution is slightly larger than the 83° contact angle found in the 
present study and in references 3 and 5. 

The U 00 curves reach a peak {c^69)max near the end of the contact arc. 
This peak was represented by a tensile stress-concentration factor In 

references 6 and 7, it was assumed that 9^ = 90°, and K^b was calculated 
as 0.81. In references 3 and 5, however, K^b was calculated as 0.92 and 0.91, 
respectively, when 0^ = 83° was used. As shown in figure 3, these results 
bracket the present K^b 0.87. In general, the present results agree 
reasonably well with the available solutions for pin-loaded holes in isotropic 
sheets. 

As an additional check on the present finite-element model, stresses were 
computed for an open hole (zero bolt modulus) in a sheet under remote uniaxial 
tension. For isotropic material properties, a gross-section stress-concentration 
factor K^g of 3.01 was calculated, compared with the familiar value of 3.00 for 
an infinite sheet. For the extreme case of 0® laminate properties, K^g = 6.45, 
which compares favorably with the infinite-sheet value of 6.43 from reference 8 . 
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width Effects 


As previously mentioned, six laminates were analyzed. The w/d values 
ranged from 2 to 20 and e/d was 10. For each case, the laminate was loaded 
by a steel pin having a snug fit. Results are first presented as stress dis- 
tributions near the hole and are briefly discussed in terms of the basic fail- 
ure modes: net-tension, shear-out, and bearing. Next, stress-concentration 

factors are plotted for a range of w/d values. Finally, these stress- 
concentration factors are compared with their corresponding infinite- laminate 
values. The quasi- isotropic laminate serves as a reference case to evaluate 
the effects of anisotropy. 

Stress distributions .- Stresses for the quasi- isotropic laminate are shown 
in figure 4. For simplicity, only results for w/d =2, 3, and oo are shown; 
other results for the full w/d range are summarized in table II. In fig- 
ure 4(a), the ctqq curve for w/d = 2 has a maximum value (^00)max that is 
more than twice that for w/d = oo. in contrast, the cXj-j- curves show only 
minor differences for the various w/d values. In fact, the maximum is 

nearly unaltered, and the contact angle 9^ varies only from about 83^ to 90^. 
Figure 4(b) shows the shear stress computed for the shear-out failure 

plane, y = R. The narrow laminate width of w/d = 2 elevated this by 

about 25 percent above the infinite-laminate case. 

The CT00, ^rr' ^xy figure 4 are associated, respectively, with 

the following failure modes: net- tension, shear-out, and bearing. Of these, 

only CJQQ shows much dependence on w/d. As a result, only the net-tension 
failure strength should be strongly influenced by laminate width. 

The 0^ laminate stresses are presented in figure 5. The ctqq curves for 
this highly orthotropic case have larger maximum values and larger gradients 
than for the quasi-isotropic case. However, (^00) max ^ only 

about twice as large as the corresponding value for w/d = oo; this w/d effect 
is comparable to that found for the quasi-isotropic laminate. The curves 

differ markedly from the cosine distribution and have a sharp peak at 6 = 0 ^, 
which is the stiff est direction for this 0^ laminate. The curves have 

almost no dependence on w/d and show a contact angle 9^ of about 82°. The 
Tj^y curves are also only slightly different for the range of w/d values. 

Results for the 90° laminate are presented in figure 6 . Although the 099 
curves have the expected peaks near 9 = 90°, they have higher peaks at 9=0°. 
Note that in this 90° laminate the fibers are tangent to the hole at the point 
of maximum CI 99 . Again, the Oj-r curves are quite similar to one another; all 
deviate from the cosine distribution. Also, they deviate significantly from the 
two previous laminates by having their maxima near 9 = 70°. The contact angles 
in figure 6 range from about 84° to 86 °. The T^^y curves in figure 6 have 
higher peaks than those found for the two previous laminates. 

The [0°/90°]g laminate results are plotted in figure 7. As expected, these 
results are approximately equal to averages of those shown for 0® and 90° lami- 
nates. The (^00) max value for w/d =2 is located at 9 = 90° and is about 
twice as large as its infinite-laminate value. The curves and the Tjjy 

curves are like those for the quasi- isotropic case in figure 4. 
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The stress distributions for the [±45]g laminate, plotted in figure 8 , are 
unlike those in previous figures. The oqq curve for w/d = has a pro- 
nounced peak at 0 = 45^, where the ~45^ ply fibers are tangent to the hole. 

For smaller w/d values, the 099 curves have a second peak near 9 = 135^ 
where the 45® ply fibers are tangent to the hole. This second peak appears to 
be strongly amplified for narrow laminates. Furthermore, narrow laminates show 
larger contact angles and larger stresses near the end of the contact arc than 
wide laminates. 

The [0®/±45^]g results shown in figure 9 have the expected trends. For 
small w/d values, (^90) max nearly three times the inf inite-* laminate value. 
This peak value occurs slightly beyond the region of contact. The and 

curves are similar to the quasi- isotropic results and are not strongly 
influenced by the laminate width. 

Stress-concentration factors .- To compare the six different laminates, the 
stress-concentration factor was calculated using each value of (<^ 00 ) max 

in the previous figures. These K^b values are listed in table II along with 
the location for each value of (^ 00 ) max contact angle 6 ^ that 

occurred for each case. The K^b values are also plotted in figure 10 for the 
range 2 i w/d ^ 10. This figure shows similar width effects for all six lami- 
nates; the larger w/d values produced smaller K^b values. As previously 
mentioned, the quasi-iso tropic case provides a reference for evaluating the 
anisotropy effects. The curve for this reference case lies near the lower 
extreme for the range of computed results. Only the [0®/±45®]g curve is lower. 
Note that the results from reference 9 for a finite-width isotropic sheet agree 
closely with the present quasi-isotropic results. As expected, the 0^ laminate 
curve is the upper bound for the range of results, and the [0®/90®]3 curve lies 
near the middle of this range. Comparisons involving the 90® laminate and the 
[ ±45®] 5 laminate are difficult because their P00)max values occurred at 
0j^ = 0® and 45®, respectively, rather than at the usual 9jj^ =« 90® for the 
laminates containing some 0 ® plies. 

The curves in figure 10 appear to have a negative slope even at w/d = 10, 
suggesting that the width effects may exist beyond w/d - 10. To investigate 
this, the Kbb values were each normalized by their corresponding infinite- 
laminate value (Ktb)oo then replotted in figure 11. Of the six curves in 

figure 11, only those for the 90® laminate and the [ ±45®] g laminate approach 
unity, indicating that only for these laminates is the effect negligible for 
w/d = 10. The other four curves converge to about 1.17. This means that 
is amplified by 17 percent even for a rather wide plate with w/d = 10. For 
comparison, if the present laminates had unloaded holes and remote uniaxial 
loads, their stress concentrations would be amplified by less than about 2 per- 
cent at w/d = 10 (ref. 10). This comparison shows that a finite width influ- 
ences the tensile stress concentration at a loaded hole much more than it does 
for an unloaded hole. In figure 11, the curves for all the laminates that con- 
tain some 0® plies lie close to the quasi- isotropic reference curve. In fact, 
for w/d ^ 3 , use of the quasi-isotropic curve to approximate the others with 
0® plies would result in errors of less than 10 percent. The other two lami- 
nates, 90® and [ ±45®] g, have curves in figure 11 that lie well below that for 
the quasi- isotropic case. However, these tw laminates have their peak values 
at 0 JJJ = 0 ^ and 45^, respectively. 
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Edge-Distance Effects 


As previously described, the six laminates were analyzed for edge-distance 
effects by using a value of w/d = 20. The e/d values ranged from 1 to 20. 

The laminates were assumed to be loaded by a steel pin fitting snugly in the 
hole. First, the results are plotted as stress distributions for each laminate. 
The stress-concentration factors for the six laminates are then compared with 
one another cuid with their corresponding infinite-laminate values. 

Stress distributions ,- The quasi- isotropic case was again used as a refer- 
ence to investigate the effects of laminate anisotropy on stresses at a loaded 
hole. Stress distributions for the quasi- isotropic case are presented in fig- 
ure 12. The curves for two finite edge distances (e/d = 1 and 2) are shown 
along with those for an infinite laminate (e/d = oo) . For the extreme e/d = 1 
case in figure 12 (a), the peak O 0 q is about twice as large as that for 
e/d = oo. As expected, this peak occurs near the end of the contact arc. But 
for small e/d values, this peak shifts slightly to the left, in contrast to 
figure 4(a), which shows a slight shift to the right for small w/d values. 
Otherwise, however, the oqq trends in figure 12(a) are quite similar to those 
in figure 4(a). The and curves in figures 12(a) and 12(b) show 

little influence of e/d, as long as e/d ^ 2 . 

Figure 13 presents stresses computed for the 0® laminate. This extremely 
orthotropic laminate is strongly influenced by finite edge distances. Even for 
e/d = 2, figure 13(a) shows that (^00) max about twice as large as for the 
reference case (e/d = oo) . For e/d = 1, it is about three times as large. 

(See table III.) This stress amplification due to finite e/d values is larger 
than that shown in the previous section for finite w/d values. This is prob- 
ably attributable to the fact that the 0 ® laminate is very stiff in the edge- 
distance direction but very flexible in the width direction. In contrast to 
the Oqq curves, the a^-^. curves show very little e/d influence; the contact 
angle 9^ is about 82^ for the full range of e/d values. 

Figure 13(b) shows that e/d has a strong influence on the values 

computed for the shear-out plane (y = R) . The x^^y curve for e/d = 1 has 
a peak that is about three times as large as the inf ini te- laminate value. The 
(Xxy)max values are similar for both the 0^ and quasi-isotropic laminates. Of 
course, these laminates have different shear strengths and therefore should 
respond differently to the same X™ stress. For example, the 0^ laminate 
has a shear strength that is only about one-fourth that for the quasi- isotropic 
case and is therefore more vulnerable to the shear-out failure mode. 

The 90® laminate stresses are plotted in figure 14. All three stresses 
show relatively little influence of e/d. These stresses were increased by only 
about 50 percent even for the extreme e/d = 1 case. Perhaps this laminate 
shows little sensitivity to the free edge because it has a very low stiffness 
in the edge-distance (x-axis) direction. 

The [0®/90®]g laminate is relatively stiff in the x-direction and, as 
expected, has a much stronger dependence on e/d than observed for the 90® 
laminate. For e/d = 1 (fig. 15), the (CJ00)max about twice as large 
as when e/d = oo. Also, (c^00)max occurs at 0 = 0® for this e/d = oo 
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reference case but occurs at 0 = 90^ for finite edge distances. Similarly, 
the peak occurs at 0 = 0 *^ for e/d = oo but shifts toward 0 = 45 ® for 

finite e/d values. 

Figures 16 and 17 present stress distributions for the [±45® ]s and the 
[0®/±45®]s laminates, respectively. As expected, the [±45®]g curves in fig- 
ure 16(a) have their peak value at 0 = 45®, where the -45® ply fibers are 
tangent to the hole. This highly stressed -45® ply is probably responsible for 
the relatively large peak shown in figure 16(b). The 099 curves for 

the [0®/+45®]g laminate have their peak at the usual 0 = 90® location. In 
fact, all the [0®/±45®]g results in figure 17 are quite similar to the 
[0®/±45®/90® ]g reference case in figure 12. Apparently, the 90® plies have 
little influence on the stresses near a loaded hole when they are used together 
with 0® and ±45® plies- 

S tress-concentration factors .- Stress-concentration factors ^or the 

six laminates are compared for a range of e/d in figure 18. These K^b values 
were calculated using (c^09)max 099 stress distributions. These K^b 

values are also presented in table III, along with the contact angle 0^ and 
the angular location for each case. Figure 18 shows that, despite the wide 

range of K^b e/d = Tf ell six curves lie within a range of 0.90 to 2.16 

for e/d =10. As in the previous analysis for w/d effects, the quasi- 
isotropic reference curve lies near the lower end of this range. Apparently, 
the quasi- isotropic laminate provides the most uniform load transfer near a 
loaded hole. Anisotropy produces more concentrated load paths and correspond- 
ingly higher stress concentrations. 

Figure 19 presents a comparison of these K^b results in a different way. 
Each Ktb value was normalized by the corresponding (Ktb)oo value for that 
laminate. Each curve in figure 19 should approach unity for sufficiently large 
e/d values, and they all appear to have this trend. For e/d = 10, each curve 
lies within the narrow range of 1.0 to 1.1. However, for e/d = 3, a practical 
value for bolted joints, the K^b ^®^ ^ quasi- isotropic laminate is about 
25 percent larger than for the infinite-laminate case; K^b values for all the 
other laminates with 0® plies are even larger. The largest increase is about 
50 percent for the 0 ® laminate. 

Because shear-out is a viable failure mode for finite e/d cases, stress- 
concentration factors were also calculated for the T™ stress along the shear- 
out planes. Figure 20 shows the stress-concentration factor Kgb = (“^xy^max/^b* 
(Also see table IV.) For e/d = 1, the Kgb results are rather tightly grouped 
in the range 0.7 to 0.9. But for larger e/d values, Kgb ranges from about 
0.2 to 0.6. As expected, the 90® laminate provides the upper bound for this 
range and the 0® laminate provides the lower bound. The quasi- isotropic curve 
lies near midrange. 

Each Kgb was normalized by the corresponding infinite-laminate value 
(Ksb)oo* These Kgb/(Kgb)oo results are shown in figure 21. As expected, for 
small e/d values, the results have a wide range. But for e/d ^ 3, all six 
curves are closely grouped and approach unity. For e/d ^ 3, the e/d effect 
is less than about 25 percent, even for the extremely orthotropic 0® laminate. 
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Stress-concentration factors for the stress were also calculated, 

but the e/d effects were found to be quite small. For e/d ^ 3, the peak 
Cj-j. stresses were less than 10 percent larger than corresponding infinite- 
laminate values. Consequently, e/d should have a negligible effect on bear- 
ing strength for e/d ^3. 


CONCLUDING REMARKS 

Stresses were calculated for finite-size orthotropic laminates that are 
loaded by a frictionless steel pin in a circular hole of the same diameter. 

These calculations were based on finite-element analyses for six laminates: 
quasi- isotropic [ 0^/145^/90^] g, 0^, 90°, [0°/90°]g, [±45°]g, and [0°/±45°]g, 
Stress-concentration factors were based on nominal bearing stress and were 
determined for wide ranges of width-to-diameter ratio w/d and edge-distance- 
to-diameter ratio e/d. For reference purposes, the infinite-laminate case 
was also analyzed for each laminate. 

When load is introduced into a laminate by a pin, the hole surface and 
pin remain in contact only over a portion of the hole boundary defined by a 
contact arc. In general, for cases with either an initial clearance fit or an 
interference fit, this contact arc varies with load. However, recent liter- 
ature shows that the present case with an initially snug fit has a contact arc 
that is independent of load. Although each combination of laminate material 
and configuration produced a different contact arc, each such arc was in the 
narrow range of 160° to 180°. Therefore, neither the wide ranges of material 
anisotropy nor laminate configuration had much influence on the contact arc. 

As expected, anisotropy had a strong influence on the contact stress distri- 
butions. These distributions deviated significantly from the well-known cosine 
distribution, especially for the 0° and 90° laminates. 

The anisotropy also had a strong influence on the magnitude and location 
of the peak tangential stress on the hole boundary. This peak was represented 
by a stress-concentration factor K^b- For the infinite-laminate case, the 
stress concentration factor ranged from 0.82 to 2.16 compared to 0.87 for the 
quasi-isotropic laminate. The stress peak developed near the end of the contact 
arc except for the 90° and [ ±45°] g laminates, which had their peaks where the 
ply directions (fibers) were tangent to the hole. 

The narrow laminates produced large K^b values for the loaded holes. Even 
for a rather wide (w/d = 10) configuration, the quasi-isotropic laminate had a 
K^b that was nearly 20 percent larger than its inf ini te- laminate reference value. 
In contrast, for an unloaded hole, the corresponding stress-concentration factor 
is only about 2 percent larger than its infinite- laminate value. In general, 
the orthotropic laminates containing 0° plies produced trends quite like that for 
the quasi- isotropic laminate, especially for w/d ^ 3. Also for w/d ^3, the 
peak bearing stresses for all six laminates were elevated by less than about 
10 percent. 

As expected, the finite-edge distance also increased the stresses near the 
loaded hole. However, this influence varied widely for the different laminates 
and different stresses. For the extreme case of e/d =1, this influence 
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increased the 0^ laminate by more than 200 percent but increased the 

90° laminate K^b by only 50 percent. The edge distances only weakly 
influenced the peak values for shear-out and bearing stress. For e/d ^ 3r 
the peak shear stresses were increased by less than about 25 percent. Simi- 
larly, for e/d i 3, the peak bearing stresses were increased by less than 
10 percent. 

In sunmary, both finite width and finite edge distance had a strong influ- 
ence on the maximum tensile stress, but had much less influence on the other 
peak stresses near the hole. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
April 1, 1981 
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APPENDIX 


APPROXIMATION FOR INFINITE- LAMINATE CASE 

To analyze the infinite-laminate case, the finite-element procedure was 
systematically applied for increasing values of w/d. The computed stress- 
concentration factor Ktb 3 loaded hole was expected to asymptotically 

approach an in finite- laminate value (Ktb)oo- This ancillary "convergence" 
study was conducted with a square, quasi-isotropic laminate and, for com- 
parison with the literature, the pin was assumed to be rigid. 

Figure 22 and table V present values for the three different load- 

ing cases analyzed in the convergence study. Case 1 , represented by the upper 
curve, is the basic problem dealt with in this study. Case 2 is included for 
comparison purposes. Note that cases 1 and 2 have the remote stress reactions 
on opposite ends of the model. Otherwise, these two cases are identical and, 
for an infinite laminate, they should have the same (^tb^ oo value because the 
remote stress would be zero for both cases. On the other hand, the difference 
between the K^b curves for cases 1 and 2 is a measure of the finite-size 
effect for a given w/d value. The curves for cases 1 and 2 converge slowly 
in figure 22 and differ by almost 20 percent at w/d = 20 . 

Because neither the case 1 curve nor the case 2 curve approached a limit 
for the present range of w/d, the third case was considered. This case could 
be generated by superimposing cases 1 and 2. Therefore, if cases 1 and 2 have 
the same limit, then the limit must also apply to case 3. Conversely, a limit 
found for case 3 would also be the desired limit for case 1. Toward this end, 
case 3 loading was used in the finite-element analysis for a range of w/d 
values. The resulting K^b curve in figure 22 has nearly converged for 
w/d = 20 . 


To check the assumption that all three curves in figure 22 converge to 
(Ktb)oo' each curve was extrapolated by using a numerical analysis. Each curve 
was approximated by the cubic polynomial 


^tb = ao + ^1 





The coefficients aQ to 33 were obtained by a least-squares fit to the 
values in table V. The polynomials were each found to have a correlation coef- 
ficient that deviated from unity by less than 4 x 10”^, indicating nearly 
perfect agreement with the tabulated K^-j^ values. For very large w/d values, 
each polynomial approaches its constant ag term. Therefore, each ag was 
an estimate for (Ktb)oo* The three ag values were 0.863, 0.864, and 0.863, 
respectively, for cases 1, 2, and 3. Therefore, the three curves converge to 
the same limit of about 0.863. 
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APPENDIX 


For convenience, the case 3 results for w/d = 20 were assumed to 
approximate the infinite- laminate results for case 1. The resulting value 
(Ktb)< 3 o® 0.869 was considered to be sufficiently close to the 0.863 extrapolated 
limit. The infinite-laminate stress distributions and (Ktb^oo values for the 
other laminates were also approximated using case 3 loading with w/d = 20. 

The (Ktb)oo values are presented in table II and, for comparison purposes, are 
repeated in table III. The shear stress-concentration factors (Ksb^oo 
infinite laminate are presented in table IV. 
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TABLE II.- TENSION STRESS-CONCENTRATION FACTORS FOR w/d RANGE WITH e/d = 10 


w/d 

Quasi- 

isotropic 

[0°/±45°/90°]s 

0° 

90° 

[0°/90°]g 



[±45° ]g 

[0°/±45°]g 

Ktb 

9m 

(9c) » 

deg 

Ktb 

(0c) » 
deg 

Ktb 

CD CD 

Ktb 

®m 

(0c) r 
deg 

Ktb 

(9c)' 

deg 

Ktb 

®m 

(9c) » 
deg 

2 

2.20 

94 

(90) 

4.37 

90 

(86) 

3.45 

0 

(86) 

3.33 

90 

(83) 

a 

2.20 

45 

(93) 

2.28 

98 

(89) 

3 

1.56 

90 

(87) 

3.29 

90 

(84) 

2.67 

0 

(84) 

2.49 

90 

(82) 

1.92 

45 

(85) 

1.53 

94 

(86) 

5 

1.22 

90 

(84) 

2.69 

90 

(82) 

2.24 

0 

(84) 

2.05 

90 

(82) 

1.71 

45 

(83) 

1.19 

90 

(83) 

10 


86 

(83) 

2.33 

90 

(82) 

2.03 

0 

(84) 

1.78 

90 

(82) 

1.57 

45 

(82) 

0.97 

90 

(83) 

20 

0.94 

m 

2.16 

m 

1.98 

B 

1.66 

90 

(82) 

1.53 

B 




°0.87 

m 

°2.00 

90 

(82) 

°1.97 

B 

c,d 

1.53 

90 

(82) 

°1.48 


°0.82 

86 

(83) 


a(a09)„ax = 2.40 at e = 120°. (See fig. 8.) 

^An approximation using w/d = 20; see case 3 loading in appendix. 
^Ktb = (Ktb)co- 

‘^('^e9)max = 1.65 at 9 = 0°. (See fig. 7.) 



























TABLE III.- TENSION STRESS-CONCENTRATION FACTORS FOR e/d RANGE WITH w/d 


Quasi- 

isotropic 

[0°/±45°/90°]g 


[0°/90O], 



[±4 

5°]g 

Ktb 

®in 

(0c)' 

deg 

2.90 

45 

(81) 


[0°/±45°], 


(0c) » 
deg 



’= 2.00 


t>1.97 




.48 


t>0.82 


An approximation using e/d = 10; see case 3 loading in appendix. 
Ktb = (Ktb)co- 

(C^ee)max = 1.65 at 0 = 0°. (See fig. 7.) 




































TABLE IV.- SHEAR STRESS-CONCENTRATION FACTORS FOR SHEAR-OUT PLANE WITH e/d RANGE AND w/d = 20 


e/d 

Quasi- 
isotropic 
■ [ 0°/±45°/90°] s 

0° 

90° 

[ 0°/90°] s 

[±450] 3 

[0°/±45O]g 

Ksb 

‘ (x/R)® 

! Ksb 

x/R 

Ksb 

! x/R 

Kgb 

x/R 

! Kgb 

x/R 

Ksb 

x/R 

1 

0.80 

0.75 

0.70 

0.80 

0.78 

0.65 

0.72 

0.80 

0.92 

0.55 

0.78 

0.60 

2 

0.55 

0.75 

0.39 

0.90 

0.64 

0.60 

0.46 

0.80 

0.66 

0.60 

0.50 

0.75 

D 

0.51 

0.75 

0.30 

1 .0 




0.80 

0.60 



in 

5 

0.47 

0.75 

0.26 

1.0 



0.37 

0.80 

0.57 

0.60 

0.43 

0.80 

10 

0.46 

0.75 

0.24 

1 .0 

0.59 

0.60 

0.35 

0.80 

0.55 

0.60 

0.42 

0.80 

b 

CO 

^^0.45 

0.75 

^0.24 

1.0 

^^0.57 

0.60 

°0.33 

0.80 

00.52 

0.65 

00.41 

0.85 


^Location of peak shear stress on shear-out plane (y = R) . 

*^An approximation using e/d = 10; see case 3 loading in appendix. 

*^Ksb = (Kgb)oo. 




































TABLE V.- STRESS-CONCENTRATION FACTORS FOR A SQUARE QUASI-ISOTROPIC 
LAMINATE LOADED BY A RIGID PIN: 0^ = 83° 


Tension 
(case 1} 


1.832 

*^(1.728) 


0m» <3eg 


Remote loading type^ 


Conpression 
(case 2) 


.576 
(0.542) 


a 

1 

j 

Tension/compression 
(case 3) 


Ktb 

<aeg 

1 


1.157 

(1.144) 


1.304 

^^( 1 . 283 ) 


.693 
(0.645) 


1.040 


0.752 

‘^(0.736) 

0.777 

^•(0.774) 

0.794 


®See figure 22 for description of loading types 
^Stress concentration at 6 = 86°. 



0.888 


0.869 




















direction 















0 90 180 

0 , deg 


(a) Stresses along hole boundary (r = R) . 



x/R 


(b) Shear stress along shear-out plane (y = R) . 
Figure 7.- Stresses in [0^/90^]g laminate as a function of w/d 



0, deg 


(a) Stresses along hole boundary (r = R) . 



x/R 


(b) Shear stress along shear-out plane (y = R) . 

Figure 8.- Stresses in [±45®]g laminate as a function of w/d; e/d 






(b) Shear stress along shear-out plane (y = R) . 


Figure 9.- Stresses in [0°/45°]s laminate as a function of w/d 














6, deg 


(a) Stresses along hole boundary. 



x/R 


(b) Shear stress along shear-out plane (y = R) 
Figure 13.- Stresses in 0® laminate as a function of e/d 





Figure 15.- Stresses in [ 0°/90] g laminate as a function of e/d 
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Figure 19.- Normalized Oqq stress-concentration factors as a function of 

e/d; w/d = 20. 









e/d 

Figure 21 Normalized stress-concentration factors for shear-out 

plane (y = R) ; w/d = 20. 
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